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The critical end-point (CEP) and critical behaviour in its vicinity, has been explored in the two 
flavour effective chiral models with and without the presence of effective Polyakov loop potential. The 
tricritical point (TCP) in the massless chiral limit has been located on the phase diagram in the /j, and 
T plane for the Polyakov loop extended Quark Meson Model (PQM) and pure Quark Meson (QM) 
model which become effective Quantum-chromodynamics (QCD) like models due to the proper 
accounting of fermionic vacuum loop contribution in the effective potential. The proximity of the 
TCP to the QCD critical end-point (CEP) has been quantified in the phase diagram. The critical 
£Sj ■ region around CEP has been obtained in the presence as well as the absence of fermionic vacuum 

t— I ' loop contribution in the effective potentials of PQM and QM models. The contours of appropriately 

normalized constant quark number susceptibility and scalar susceptibility have been plotted around 
CEP in different model scenarios. These contours determine the shape of critical region and facilitate 
^jqJ comparisons in different models such that the influence of fermionic vacuum term and Polyakov loop 

potential on the critical behavior around CEP can be ascertained in qualitative as well as quantitative 
, terms. Critical exponents resulting from the divergence of quark number susceptibility at the CEP, 

■ have been calulated and compared with in different model scenarios. The possible influence of TCP 
' on the critical behavior around CEP, has also been discussed. The temperature variation of a and 

y—( , 7r meson masses at fi = 0, fi = iicep and /i > ficEP has been shown and compared with in different 

model scenarios and the emerging mass degeneration trend in the a and tt meson mass variations 
i i has been inferred as the chiral symmetry restoration takes place at higher temperatures. 

; 

P-l PACS numbers: 12. 38. Aw, 11.30.Rd, 12.38.Lg, ll.10.Wx 

I. INTRODUCTION 

> " 

" Under the extreme conditions of high temperature and/or density, normal hadronic matter undergoes a phase 
, transition, where the individual hadrons dissolve into their quark and gluon constituents and produce a collective form 
of matter known as the Quark Gluon Plasma(QGP)[l|-|j|. Study of the different aspects of this phase transition, is 
C^l \ a tough and challenging task because Quantum Chromodynamics(QCD) which is the theory of strong interaction, 

■ becomes nonperturbative in the low energy limit. However the QCD vacuum reveals itself through the process of 
| spontaneous chiral symmetry breaking and phenomenon of color confinement. 

CN ■ The QCD Lagrangian is known to have the global SUL+n(Nf) x SUL-n(Nf) symmetry for Nf flavours of massless 
quarks. The formation of a chiral condensate in the low energy hadronic vacuum of QCD, leads to the spontaneous 
J> . breaking of the axial (A=R-L) part of this symmetry known as the chiral symmetry and one gets (Nj — 1) massless 
Goldstone bosons according to the Goldstone's theorem. Since quarks are not massless in real life, chiral symmetry of 
the QCD lagrangian gets explicitly broken and massless modes become pseudo- Goldstone bosons after acquiring mass. 
Nevertheless, the observed lightness of pions in nature suggests that we have an approximate chiral symmetry for 
QCD with two falvours of light u and d quarks. In the opposite limit of infinitely heavy quarks, QCD becomes a pure 
SU (N c ) gauge theory which remains invariant under the global Z(N C ) center symmetry of the color gauge group. The 
Center symmetry which is a symmetry of hadronic vacuum, gets spontaneously broken in the high temperature/density 
regime of QGP. The expectation value of the Wilson line (Polyakov loop) is related to the free energy of a static color 
charge. It vanishes in the confining phase as the quark has infinite free energy and becomes finite in the deconfined 
phase. Hence the Polyakov loop serves as the order parameter of the confinement-deconfinement phase transition [f|. 
Even though the center symmetry is always broken with the inclusion of dynamical quarks in the system, one can 
regard the Polyakov loop as an approximate order parameter because it is a good indicator of a rapid crossover in the 
confinement-deconfinement transition Iq, Mj . 

Lattice QCD simulations (see e.g. |8l4l8j) give us important information and insights regarding various aspects of 
the QGP transition, like the restoration of chiral symmetry in QCD, order of the confinement-deconfinement phase 
transition, richness of the QCD phase structure and phase diagram mapping. Since lattice calculations are technically 
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involved and various issues are not conclusively settled within the lattice community, one resorts to the calculations 
within the ambit of phenomenological models developed in terms of effective degrees of freedom. These models serve 
to complement the lattice simulations and give much needed insight about the regions of phase diagram inaccessible 
to lattice simulations. 

Construction and mapping of the phase diagram in the quark chemical potential and temperature plane is the prime 
challenge before the experimental as well as theoretical QGP community. On the temperature axis, the chiral transition 
at zero quark chemical potential with almost physical quark masses, has been well established to be a crossover in 
recent lattice QCD simulations [HI, Ojj]. Effective chiral model studies[2(| predict first order phase transition at lower 
tempertaures on the chemical potential axis. Thus the existence of a critical end point(CEP) has been suggested 
in the phase diagram based on model studies [21H23T ] together with the inputs from lattice simulations [Tol - fT^ . The 
first order transition line starting from the lowest temperarute on the chemical potential axis, terminates at the CEP 
which is a genuine singularity of the QCD free energy. Here the phase transition turns second order and its criticality 
belongs to the three dimensional Ising universality class ■ The precise location of the CEP is highly sensitive to 

the value of the strange quark mass. Lattice QCD predictions at non zero chemical potential are much more difficult 
due to the QCD action becoming complex on account of the fermion sign problem j9| . There is evidence for a CEP 
at finite /it[lol. Illj from a Taylor expansion of QCD pressure around /j, — 0, however in another lattice study, finite 
chemical potential extrapolations provide some limitations and can rule out the existence of a CEP for small fi/T 
ratios[28]. In the chiral limit of zero up and down quark masses, the chiral phase transition is of second order at 
zero n and the static critical behavior is expected to fall in the universality class of the 0(4) spin model in three 
dimensions 20] . Thus the existance of CEP for real life two flavor QCD implies that two flavor massless QCD has a 
tricritical point(TCP) at which the second order 0(4) line of critical points ends. 

Experimental signatures encoding the singular behavior of thermodynamic quantities in the vicinity of critical point 
have already been suggested 29] . These are related to chemical potential and temperature fluctuations in event- by- 
event fluctuations of various particle multiplicities [30j . In the centre of mass energy scans, an increase and then a 
decrease in the number fluctuations of pions and protons should be observed as one crosses the critical point. If 
the signals are not washed out due to the expansion of the colliding system, the critical point might be located in 
the phase diagram by the observation of nonmonotonic behavior of number fluctuations in its vicinity [3 1[. Recently 
"beam energy scan" program dedicated to the search of critical point has been started at the Relativistic Heavy 
Ion Collider (RHIC, Brookhaven National Laboratory) experiments [32j . The Compressed Baryonic Matter (CBM) 
experiment (GSI-Darmstadt) at the facility for Antiproton and Ion Research (FAIR) and the Nuclotron-Based Ion 
Collider facility (NICA) at the Joint Institute for Nuclear Research (JINR), will also be looking for the signatures of 
critical end point. Characteristic signatures of the conjectured CEP for experiments have been discussed in refs [33i - [35| . 

Recently, effective chiral models like the linear sigma models(LSM) [36M42l | .the quark-meson (QM) modcls(scc 
e.g.[i^-[53), Nambu-Jona-Lasinio (NJL) type models [H, [54l-l57|. were extended to combine the features of 
confinement-deconfinement transition together with that of chiral symmetry breaking-restoring phase transition. Chi- 
ral order parameter and Polyakov loop order parameter got simultaneously coupled to the quark degrees of freedom 
in these models. Thus Polyakov loop augmented PNJL models [58l - [74| ,PLSM models and PQM models [75-83:] have 
facilitated the investigation of the full QCD thermodynamics and phase structure at zero and finite quark chemical 
potential and it has been shown that bulk thermodynamics of the effective models agrees well with the lattice QCD 
data. The issue of location of CEP in phase diagramn tog ether with the extent of criticality around it, is also being 
actively pursued in a variety of effective model studies [2J, |25j [84j-|90j . The critical region around CEP is not point- 
like but has a much richer structure. The estimation of the size of critical region is especially important for future 
experimental seraches of CEP in heavy-ion collision experiments. 

In the no-sea mean-field approximations, an ultraviolet divergent part of the fermionic vacuum loop contribution to 
the grand potential got frequently neglected till recently in the QM/PQM model calculations [43l 1441. |49|. l5ll. l56j . Due 
to this, the phase transition on the temperature axis at fi = for two flavour QM model becomes first order in the 
chiral limit of massless quarks and one does not find TCP on the phase diagram. Recently, Skokov et al. in Ref. 



addressed this issue by incorporating appropriately renormalized fermionic vacuum fluctuations in the thermodynamic 
potential of the QM model at zero chemical potential which becomes an effective QCD-like model because now it 
can reproduce the second order chiral phase transition at /i = as expected from the universality arguments [20j for 
the two massless flavours of QCD. The fermionic vacuum correction and its influence has also been investigated in 
earlier works[;9ll494|. In a recent work [95], we generalized the proper accounting of renormalized fermionic vacuum 
fluctuation in the two flavour PQM model to the non-zero chemical potentials and found that the position of CEP 
shifts to a significantly higher chemical potential in the /i and T plane of the phase diagram, due to the influence of 
fermionic vacuum term in our PQMVT (PQM model with vacuum term) model calculations. Very recently, Schaefer 
ct. al.[96] worked out the size of critical region around CEP in a three flavour (2+1) PQM model where cut off 
independent renormalization of fermionic vacuum fluctuation has been considered. They calculated critical exponents 
and higher order non-gaussian moments to identify the fluctuations in particle multiplicities. Since the criticality 
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around CEP is influenced by the presence of strange quark, it is important to have a two flavor calculation in the 
same model in order to faclitate the comparion with the corresponding size of critical region and nature of criticality 
obtained in 2+1 flavour QM/PQM model studies. 

In this paper, we will calculate the phase diagram in the massless chiral limit and locate the tricritical point (TCP) 
in the fj, and T plane for the PQMVT and QMVT (QM model with vacuum term) models which have become QCD- 
like in the presence of fermionic vacuum term and yield the second order transition at \i = on the temperature 
axis. Further, we will be investigating the size and extent of critical region around the CEP in phase diagram 
calculated in the two flavour QM /PQM models with and without the effect of fermionic vacuum fluctuations in the 
grand potential. We will be plotting the contours of appropriately normalized constant quark number susceptibility 
and scalar susceptibility around CEP in different model scenarios. In order to investigate the qualitative as well 
as quantitative effect of fermionic vacuum term and Polyakov loop potential, on the critical behavior around CEP, 
we will compare the shape of these contours as obtained in different model calculations. Further, we compute and 
compare the critical exponents resulting from the divergence of quark number susceptibility at the CEP in different 
model scenarios. The possible influence of TCP on the critical behavior around CEP, will also be discussed. Finally, 
we plot the temperature variation of a and 7r meson masses at /i = 0, fi = [icep and /i > [icep in different model 
scenarios and compare the emerging mass degeneration trend in the cr and 7r meson mass variations as the chiral 
symmetry gets restored at higher temperatures. 

In the presentation of this paper, we recapitulate the formulation of the two quark flavour PQM model in Secim 
The thermodynamic grand potential and the choice of the Polyakov loop potential has been discussed in subsection 
III Al In the subsection III Bl we give a brief description of the appropriate renormalization of fermionic vacuum loop 
contribution and explain how the new model parameters are obtained in vacuum when renormalized vacuum term is 
added to the effective potential. The section Mil explores the proximity of QCD tricritical point to the critical end- 
point and the detail structure of the phase diagram for the QMVT and PQMVT models where the effect of fermionic 
vacuum term has been taken care of in the QM and PQM models. The structure of the phase diagram for QM and 
PQM model and the location of critical end point has also been presented to facilitate the comparison. The subsection 
1111 Al investigates the extent of criticality around CEP where contours of constant baryon number susceptibility ratios 
and constant scalar susceptibility ratios, have been presented in the p, and T plane and comparison in all the four 
models QM,PQM, QMVT and PQMVT, have been made. The critical exponents for the criticality around CEP in 
all the four models QM,PQM, QMVT and PQMVT, have been discussed in the subsection HITBl Subsection HTTP 
presents the temperature variation of a and n meson masses at \i = 0, fi = (J-cep and fi > \icep- Here we also present 
a detail comparison of the emerging mass degeneration trends in the a and 7r meson mass variations in different model 
scenarios as the chiral symmetry restoration takes place at higher temperatures. In the end Sec. II VI presents summary 
together with the conclusion. The first and second partial derivatives of U\ og and with respect to temperature 
and chemical potential has been evaluated in appendix A of Ref. [95[ • 



II. MODEL FORMULATION 



We will be working in the two flavor quark meson linear sigma model which has been combined with the Polyakov 
loop potential [75(In this model, quarks coming in two flavor are coupled to the SUl(2) x SUr{2) symmetric four 
mesonic fields a and if together with spatially constant temporal gauge field represented by Polyakov loop potential. 
Polyakov loop field $(af) is defined as the thermal expectation value of color trace of Wilson loop in temporal direction 
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-Tr c L, 



*' = i-Tr ci t 



where L(x) is a matrix in the fundamental representation of the SU C (3) color gauge group. 



L(x) = Vexp 



drAo (x, t) 



(1) 



(2) 



Here V is path ordering, A a is the temporal component of Euclidean vector field and f3 — T 1 [j| . 

The model Lagrangian is written in terms of quarks, mesons, couplings and Polyakov loop potential U ($, $*, T). 



Cpqm = Cqm -U(<S>,$>\T) 
where the Lagrangian in quark meson linear sigma model 

£qm = <Tf [il^D^ - g(a + i-f b f ■ n)] q f + C m 



(3) 



(4) 
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The coupling of quarks with the uniform temporal background gauge field is effected by the following replacement 
= df, - iA^ and = S^qAq (Polyakov gauge), where A^ = g s A*\ a /2. g s is the SU C (3) gauge coupling. A a are 
Gell-Mann matrices in the color space, a runs from 1 • • • 8. qj — (u, d) T denotes the quarks coming in two flavors and 
three colors, g is the flavor blind Yukawa coupling that couples the two flavor of quarks with four mesons; one scalar 
(a, J p — + ) and three pseudo scalars (7?, J p = — ). 

The quarks have no intrinsic mass but become massive after spontaneous chiral symmetry breaking because of non 
vanishing vacuum expectation value of the chiral condensate. The mesonic part of the Lagrangian has the following 
form 

C m = l^a) 2 + ~(9 m tt) 2 -U(a,7f) (5) 

The pure mesonic potential is given by the expression 

U{a,ir) = ^(a 2 + Tr 2 -v 2 ) 2 -ha, (6) 

Here A is quartic coupling of the mesonic fields, v is the vacuum expectation value of scalar field when chiral symmetry 
is explicitly broken and h =f 7r m 2 . 



A. Polyakov loop potential and thermodynamic grand potential 

The effective potential U ($, $*,T) is constructed such that it reproduces thermodynamics of pure glue theory on 
the lattice for temperatures upto about twice the deconfinement phase transition temperature. In this work, we are 
using logarithmic form of Polyakov loop effective potential [58]. The results produced by this potential is known to 
be fitted well to the lattice results. This potential is given by the following expression 



^i og ($,$*,r) 

ji4 



a(T) 
2 

+4($* 3 



$*$ + 6(T) ln[I-6$*$ 
$ 3 ) - 3($*$) 2 ] 



(7) 



where the temperature dependent coefficients are as follow 



a{T) = a + ai 



a 2 



The parameters of Eq. (O arc 



ao 
a-2 



3.51 , 
15.2 , 



a 1 
b 3 



b{T) 



-2.47 
-1.75 



The critical temperature for deconfinement phase transition To = 270 MeV is fixed for pure gauge Yang Mills 
theory. In the presence of dynamical quarks To is directly linked to the mass-scale Apnri, the parameter which has 
a flavor and chemical potential dependence in full dynamical QCD and To — > To(Nf,fx) [75], |83|. For our numerical 
calculations in this paper, we have taken a fixed T = 208 for two flavours of quarks. 

In the mean-field approximation, the thermodynamic grand potential for the PQM model is given as (75| 



tt MF (T, ft a, $, $*) = U{T- $, $*) + U(a) + fi, ff (T, W a, $, $*). (8) 

Here, we have written the vacuum expectation values (a) = a and (7?) = 
The quark/ ant iquark contribution in the presence of Polyakov loop reads 



Qvac 



at 



(9) 



The first term of the Eq. @ denotes the fermion vacuum contribution, regularized by the ultraviolet cutoff A. In 
the second term g+ and g~ have been defined after taking trace over color space. 
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9t = 



1 + 3$e-^ + / T + Z<S>*e- 2E H T + e^ E t l T 



(10) 



1 + 3$*e-^/ T + 3<$>e~ 2E ^ T + e - 3E ^ T 



Here we use the notation = E q =p /i and i? 9 is the single particle energy of quark/antiquark. 



(11) 



E q = Jp 2 + m q 2 



(12) 



where the constituent quark mass m q — ga is a function of chiral condensate. In vacuum er(0, 0) = cto = f-n = Qi.QMeV 



B. The renormalized vacuum term and model parameters 



The fermion vacuum loop contribution can be obtained by appropriately renormalizing the first term of Eq. ^ 
using the dimensional regularization scheme, as done in Ref.[82j. A brief description of essential steps is given below. 
Fermion vacuum term is just the one-loop zero temperature effective potential at lowest order [97j 



Qvac 



~ 2N ^ I (0^ 



= -2N f N c J Hp 2 + El) + K, 



(13) 



the infinite constant K is independent of the fermion mass, hence it is dropped. 

The dimensional regularization of Eq. (TIB")) near three dimensions, d ~ 3 — 2e leads to the potential up to zeroth 
order in e as given by 



^ var N c N f 4 



'II 'I 



167T 



2 9 



- 3 + ^ + 41 °(2^m) 



(14) 



here M denotes the arbitrary renormalization scale. 

The addition of a counter term 8C in the Lagrangian of the QM or PQM model 



SC = To^ V {7 - ^[-3 + 27. - 41n (2^)] } , 



(15) 



gives the renormalized fermion vacuum loop contribution as 



Now the first term of Eq. ^ which is vacuum contribution will be replaced by the appropriately renormalized 
fermion vacuum loop contribution as given in Eq. (|16p . 

The relevant part of the effective potential in Eq. ([8} which will fix the value of the parameters A and v in the 
vacuum at T — and (J, = is the purely a dependent mesonic potential U(a) plus the renormalized vacuum term 
given by Eq. (fTl)]) . 

0(0-) = n g f + U{a) = -^ 3 V In (g) - ^a 2 + V ha, (17) 
The first derivative of Q(a) with respect to a at a = f n in the vacuum is put to zero 

an M F(o,Q;<7,$,$*) on(a) 

da = ~da- =Q (18) 

The second derivative of VL(a) with respect to a at er = f v in the vacuum gives the mass of a 



6 



ad 2 



9a 2 



(19) 



Solving the equations (fl"8|) and (fT9|) . we obtain 



jyyVj 4 

A = A * + ^r* 



3 + 41n 



M 



(20) 



and 



N r N 



Xv 2 = (Xv 2 ) s + ^g*f 2 
where A s and (Xv 2 ) s are the values of the parameters in the pure sigma model 



(21) 



Xs = 



2/ 2 



(22) 



(Xv 2 



m 2 — 3m 2 



(23) 



It is evident from the equations (j2"0")) and (|2"TT) that the value of the parameters A and v 2 have a logarithmic 
dependence on the arbitrary renormalization scale M. However, when we put the value of A and Ai; 2 in Eq. (|17[) . the 
M dependence cancels out neatly after the rearrangement of terms. Finally we obtain 



8^ 2 y \fj 2 4 



—a — ha, 



(24) 



Here, we define A r and X r v 2 as the values of the parameters after proper accounting of the renormalized fermion 
vacuum contribution. 



A r — A s 



3N r N 



8tt 2 



and 



; \ .2 , , N ? N f „± r2 

4tt 2 3 /w 



A r u 2 = (Xv 2 ) s + 



(25) 



(26) 



Now the thermodynamic grand potential for the PQM model in the presence of appropriately renormalized fermionic 
vacuum contribution (PQMVT model) will be written as 



n MF (T, (j,; a, $, $*) = U{T- $, $*) + 0(<r) + fi* (T, m a, $, $ 



(27) 



Thus in the PQMVT model, One can get the chiral condensate a, and the Polyakov loop expectation values $, $* 
by searching the global minima of the grand potential in Eq. (|27[) for a given value of temperature T and chemical 
potential /i 



MF 



MP 



MF 



da 



9$ 



<9$* 



= , 



(28) 



We will take the values m, = 138 MeV, m a = 500 MeV, and /„• = 93 MeV in our numerical computation. The 
constituent quark mass in vacuum mj = 310 MeV fixes the value of Yukawa coupling g = 3.3. 
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FIG. 1: (a) For calculations with experimental pion mass, solid lines representing the first order chiral phase transition merge 
with the dotted lines (blue in color) for the chiral crossover at the CEP which is denoted by the filled circle. The thick solid 
lines around CEP are the contours of constant normalized quark number susceptibility R q =2. For calculations with zero pion 
masss, dash lines represent the first order phase transition in the chiral limit of QMVT and PQMVT models while dash dot 
lines represent the second order transition and the filled triangle is the location of TCP where these two lines merge into each 
other. Lower part of the Fig. shows the QMVT model results while upper part shows the PQMVT results, (b) Lower part of 
the Fig. shows the QM model results while upper part shows the PQM results. The line types represent the same thing as in 
Fig. a for calculations with experimental pion mass. 



III. THE PROXIMITY OF THE TCP TO THE CEP AND THE PHASE STRUCTURE 



The presence of CEP in the \x and T plane of the phase diagram for the real life two flavor QCD, implies the 
existence of a tricritical point (TCP) for the corresponding two massless quark flavor QCD because the chiral phase 
transition on the temperature axis turns second order at zero fi in the chiral limit. The second order 0(4) line of 
critical points, starting from /i=0 and finite T, ends at the TCP in the /x and T plane where it happens to meet the 
first order transition line originating from the chemical potential axis at the lowest temperature. The two flavour 
QMVT and PQMVT models, where the effect of fermionic vacuum fluctuation has been incorporated in the effective 
potential of QM and PQM models, are effective QCD-like models. Hence one must find a TCP in the fi and T 
plane of the phase diagram computed in the chiral limit of zero pion mass in these models. The present work starts 
with the computation of phase diagram and the location of the CEP in the fi and T plane of all the four models 
QMVT, PQMVT, QM and PQM for the real life explicit chiral symmetry breaking with the experimental value of pion 
mass. Next, we locate the TCP in our calculation and quantify its proximity to the CEP in the phase diagram. The 
presence or absence of TCP in the phase diagram of a model calculation and its distance from CEP, influences the 
nature of critical fluctuations around CEP. 



The results for QMVT and PQMVT model calculations with real life pion mass, have been presented in Fig 1(a) 
while Fig 1(b) | presents the corresponding results for the QM and PQM model calculations. The locations of the TCP 
in the fj, and T plane of the phase diagrams computed with zero pion mass in the QMVT and PQMVT models, have 
also been shown in Fig |l(a) The TCP does not exist in the phase diagram of QM and PQM models in the chiral 
limit of zero pion mass because the phase transition, on the temperature axis at [i = 0, has been found to be of first 
order. For calculations with experimental pion mass, solid lines representing the first order chiral phase transition in 
Fig[T] merge with the dotted lines (blue in color) for the chiral crossover at the CEP (denoted by filled circle) . The 
±5 MeV error bars (in a range // = 100 to \i = 160 MeV) on the dotted line in the upper part of Fig 1(a) signify 
the ambiguity of pseudo-critical temperature determination for the chiral crossover transition in the PQMVT model 
(see[95| for details) calculations. The thick solid lines around CEP are the contours of constant ratio (R q =2) of quark 
number susceptibility obtained in a model calculation to the value of quark number susceptibility for a free quark 
gas. Since quark number susceptibility diverges at the CEP, such contours signify the extent of critical fluctuations 
around CEP. The CEP in the QMV T mo del is located at h C ep=299.35 MeV and T C £;p=32.24 MeV as shown by the 
filled circle in the lower part of Fig 1(a) It shifts to the higher value on the temperature axis at Tc7pp=83.0 MeV 



and /ic7Bp=295.217 MeV in PQMVT model due to the influence of Polyakov loop potential. 

The filled circle in the lower part of Fig |l(b)| locates the CEP in QM model at T C ep=102.09 MeV and fj, C EP=151-7 
MeV and again in the influence of Polyakov loop potential, the CEP in PQM model shifts considerably towards the 
temperature axis at Tcpp=166.88 MeV and ^cep— 81-02 MeV in the upper part of the Fig 1(b) If we compare 
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the location of CEP in QM and PQM models as shown in Fig]l(b)| to the location of CEP in QMVT and PQMVT 



models in Fig 1(a) we find a considerably significant shift of CEP to large chemical potential and small temperature 



values for QMVT and PQMVT models due to the robust influence of fermionic vacuum term inclusion in the effective 



potential. Another important thing worth noticing for the PQMVT model phase diagram in the upper part of Fig 1(a) 



is the highest value of temperature for the chiral crossover transition occurring on the temperature axis at fi = 0. 
The combined effect of Polyakov loop potential and fermionic vacuum term is responsible for this. If one compares 



the phase diagram in the lower part of Fig 1(b) for QM model with the phase diagram in the lower part of Fig 1(a) 



for QMVT model, one immediately notices that the the chiral crossover transition at fx = occurs at a higher 
temperature value only due to the influence of fermionic vacuum fluctuation. These results are the extension of our 
recently reported work[95| and facilitate the details of model comparison for the t wo q uark flavour case. Further these 
results are also in qualitative agreement with the recent results of Schaefer et. al.[96| for the 2+1 flavour case. 

For calculations with zero pion masss, dash lines represent the first order phase transition in the chiral limit of 
QMVT and PQMVT models in Fig 1(a) while dash dot lines represent the second order transition and the filled 
triangle is the location of TCP where these two lines merge into each other. In the upper part of the Fig 1(a) the 
filled triangle locates the presence of tricritical point (TCP) at 7W37.09 MeV and ^=240.14 MeV for PQMVT model 
calculation. In order to quantify the proximity of TCP to the CEP, we have plotted the constant normalized quark- 
number susceptibility (R q =2) contour around CEP. It is seen on the phase diagram that the range and extension of 
this contour is quite large in both the directions;chemical potential as well as the temperature. The second cumulant of 
the net quark number fluctuations on this contour is double to that of the free quark gas value and such enhancements 
are the signatures of CEP for the heavy-ion collision experiments. The TCP location is quite well inside this contour 
on the phase diagram. It means that the shape of the critical region and nature of criticality around CEP, gets 
influenced by the presence of TCP in the corresponding chiral limit. In a recent NJL/PNJL model calculation by 
Costa et. al.poj. the CEP lies closer to the chemical potential axis but the TCP gets located on the periphery of 
R q =2 contour around CEP. In the QMVT model calculation, the tricritical point (TCP) is found at T 4 =69.06 MeV 
and fi t =263.0 MeV as denoted by filled triangle in the lower part of the Fig 1(a) Here also the TCP lies quite well 
inside the R q —2 contour on the phase diagram. 



Susceptibility Contours and Criticality 



In order to locate the CEP in heavy-ion collision experiments, one requires the quantification of criticality around 
CEP. The crossover transition is marked by a peak in the quark number susceptibility which becomes sharper and 
higher as one approaches the CEP in the phase diagram from the crossover side and finally the peak diverges at CEP. 
Hence the quark number susceptibilities and scalar susceptibilities will be significantly enhanced in a region around 
the CEP in the /i and T plane in comparison to their respective values for the free quark gas. Thus the contour 
regions of properly normalized constant quark number susceptibilities and scalar susceptibilities, can be taken as the 
measure of criticality around CEP. The ratio of quark-number susceptibility \q normalized to the free susceptibility 



free 



Ra 



Xq 
v frcc 



The expression of quark number susceptibility is obtained as 

Xq 



d 2 n 



MF 



m q ^0 O 

u q = 2N c N f = 12 



3 M 2 



free 
Aq 



(29) 



(30) 

(31) 
(32) 



The first and second partial derivatives of er, $ and <E>* fields with respect to chemical potential contribute in the 
double derivatives of f2(cr), U\ a% and with respect to chemical potential as given in the appendix A of Ref . [95l| . 
Contours with three different values for the ratios R q , have been plotted in Fig|2] in the [i and T plane relative 



to the CEP. If we compare the contours in Fig 2(a) depicting the PQM model results to the contours in Fig 2(b) 
showing the pure QM model results, we conclude that the presence of Polyakov loop potential, compresses the critical 
region particularly in the T direction similar to findings of Schaefer et. al.[96| in their three flavour calculation. The 
compression of critical region in the T direction is much more pronounced in our two quark flavour calculation as can 
be seen in the spread of R q = 2 contour on the temperature axis only in a small range of ±2.5 MeV near Tcep- The 
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FIG. 2: (a) The contours of three different values for the constant ratios R q = 2, 3 and 5 of quark number susceptibility to the 
quark susceptibility for the free quark gas, are plotted in the PQM model calculations. (b) Similar contours plotted in the QM 
model calculations. 
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FIG. 3: (a) The contours of three different values for the constant ratios R q = 2, 3 and 5 of quark number susceptibility to the 
quark susceptibility for the free quark gas, are plotted in the PQMVT model calculations. (b) Similar contours as plotted in the 
QMVT model calculations. 



modification in the /i-direction is quite moderate compared to the effect in the T direction. Since the chiral crossover 
transition becomes faster and sharper due to the presence Polyakov loop contribution in the effective potential, the 
critical region in the T direction gets significantly compressed. 

The size of the critical region is significantly influenced by the incorporation of fermionic vacuum fluctuations in 
the effective potential as shown in FigEl In the phase diagram, the size of critical region is increased in a direction 
perpendicular to the crossover line due to the influence of the fermionic fluctuations. This effect is less pronounced in 
Fig 3(a) because of the compression of critical region width due to the presence of Polyakov loop potential contribution 
in PQMVT model while the QMVT model results of Fig |3 (b) | obtained in the absence of Polyakov loop, show a robust 
increase in the width of the critical region. However, the extent and size of critical region in the PQMVT model 
in Fig 3(a)| is noticeably larger in both the directions \x as well as T compared to that of QMVT model results as 
shown in Fig |3(b)| In the presence of fermionic vacuum term, CEP gets located at larger chemical potentials in 
QMVT/PQMVT models. Since the quark determinant gets modified mostly at moderate chemical potentials by 
the presence of Polyakov loop potential and further in its influence, the PQMVT model CEP shifts to a higher 
critical temperature [cf. also Figfl] when compared to the CEP in QMVT model, we obtain an enhancement of the 
critical region in PQMVT model. Further the chiral crossover transition becomes much smoother because the phase 
transitions in general get washed out in the influence of fluctuations. This leads to a critical region which is broader 
in perpendicular direction to the extended first-order transition line. The influence of the Polyakov loop potential 
becomes insignificant for smaller temperatures and larger chemical potentials, hence the size of the critical region for 
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FIG. 4: (a) The contours of different values for the constant ratios Rs = 10 and 15 of T and /x dependent scalar susceptibility 
to the scalar susceptibility at T=0 and /i=0 MeV, are plotted in the PQM model calculations. (b) Similar contours for the 
constant ratios Rs = 10, 15 and 25 are plotted in the QM model calculations. 



(/i — Hcep) > becomes comparable in both the models QMVT and PQMVT. In FigEl filled circles are the position 
of CEP and the filled triangles, show the location of TCP, we observe that the TCP gets located outside R q = 3 and 
quite well inside the R q = 2 contour. 

If we compare our two quark flavour results with the 2+1 flavour calculations in renormalized PQM/QM models in 
Ref . [96| . we notice that in the absence of strange quarks, the effect of fermionic vacuum term, leads to an enhanced 
critical region in both the directions T as well as /i and the size of contours is larger in our two flavour calculation. 
We point out that the 2+1 flavour calculation in Ref.|96j was done with m a — 400 MeV and To = 270 MeV while 
in our two flavour calculation ra a — 500 MeV and To = 208 MeV. In general higher value of m a pushes the CEP to 
higher chemical potential. In our two quark flavour calculation, the CEP is at {Tcep, ^cep)—(83. 0,295. 217) MeV 
and (32.24,299.35) MeV respectively in PQMVT and QMVT model calculations while the CEP in the corresponding 
the 2+1 flavour model calculation of Ref.@ is at {T CE p, /xcep)=(90.0,283.0) MeV and (32,286) MeV. 

The zero-momentum projection of the scalar propagator, encodes all fluctuations of the order parameter and it 
corresponds to the scalar susceptibility \ a . The relation of the scalar susceptibility to the order parameter is obtained 

as m (USUI. 



The most rapid change of the chiral order parameter, should be coincident with the maximum in the temperature 
or quark chemical potential variation of Xa- The relation of scalar susceptibility to the sigma mass via Xu ~ 171 a 2 can 
be easily verified. The normalized scalar susceptibility is written as[49| 

*Cr.,)-^$ (34) 

In Fig. 21 the contours have been plotted for three values of fixed ratios R s around the CEP in the PQM and 



QM models. The R s = 10 contour in Fig. 4(a) is compressed in the T direction and its extension in /x direction is 
also reduced in comparison to the pure QM model contours in Fig. |4(b)| This is due to the quite fast and rapid 
temperature or chemical potential variation of a meson mass m a on account of faster and sharper change of order 
parameter for chiral crossover in the presence of Polyakov loop potential in the calculations. We do not find contour 



for R s = 25 in Fig. 4(a) because the minimum value of a meson mass does not fall below 100 MeV, though the value of 
m a falls very rapidly and sharply from 500 MeV to 128 MeV giving rise to a very thin and small contour region even 
for R s = 15. In the QM model calculations, we get all the contour regions for R s = 10, 15 and 25 with well defined 
size because the m a variation is smoother and slower in comparison to the corresponding PQM model results and 
further the minimum in the m a variation approaches almost zero value in QM model. The chiral crossover transition 
on the temperature axis at /i = MeV in the QM and PQM models, is quite sharp and fast because it emerges from 
the background of first order chiral transition at fj, = MeV in the corresponding chiral limit of zero pion mass and 
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FIG. 5: (a) The contours of three different values for the constant ratios Rs = 10, 15 and 25 of T and fi dependent scalar 
susceptibility to the scalar susceptibility at T=0 and ^i=0 MeV,are plotted in the PQMVT model calculations. (b)Similar 
contours for the constant ratios Rs = 10, 15 and 25 are plotted in the QMVT model calculations. 



we do not find the existence of TCP in QM/PQM models. As a consequence, we do not find the closure of R s = 10 
contour on the temperature axis at /i = MeV. 

We obtain quite well defined and closed contour regions for R s = 10, 15 and 25 in Fig. [S] which again become 
broader in the direction perpendicular to the crossover line due to the presence of fermionic vacuum fluctuations 
in QMVT and PQMVT model calculations. For scalar susceptibility also, the critical region gets elongated in the 
phase diagram and Xa is enhanced in the direction parallel to the first-order transition line. Here also the presence 
of Polyakov loop potential in the PQMVT model, leads to the compression in the width of critical region around 
CEP as shown in Fig 5(a) The fermionic vacuum fluctuations, make the chiral crossover transition very smooth 



while the Polyakov loop potential makes it sharper and faster and these opposite effects give a typical shape to the 
quark number susceptibility contours in Fig 3(a) in the PQMVT model. Similar effects can be seen in the scalar 



susceptibility contours in Fig 5(a) In the influence of fermionic vacuum fluctuations only, the \a contours in Fig 5(b) 



in the pure QMVT model are broader and rounded. 

For the detail understanding and analysis of the criticality around the CEP, we will be studying the critical exponents 
of the susceptibilities at the critical point in the next section. 



B. Critical Exponents 



The crossover transition is marked by a peak in the quark number susceptibility which diverges as one approaches 
the CEP from the crossover side in the phase diagram. This divergence is governed by a power law within the 
critical region. The corresponding critical exponents depend on the route through which the singularity (CEP) is 
approached in the /i and T plane 98] . This path dependence decides the shape of the critical region. In the mean- field 
approximation, the quark number susceptibility scales with an exponent j q = 1 for a path asymptotically parallel to 
the first-order transition line and for any other path which is not parallel to the first-order line, the divergence scales 
with the exponent e = 2/3. This larger critical exponent (7^ > e) is one reason for the elongation of the critical region 
in a direction parallel to the first-order line as already pointed out in [24. 1491] . 

In order to further investigate the nature of criticality in two flavour calculations, we have numerically evaluated 
the critical exponents of the quark number susceptibility \q m QM,PQM,QMVT and PQMVT models. In these 
investigations, the critical hcep at fixed critical temperature Tqep is approached from the lower as well as higher 
fi sides in a path parallel to the /i-axis in the (T, /i)-plane. The calculation of the critical exponents, has been done 
with the following linear logarithmic fit formula: 

logXg = -mlog|/i- ^cep\ +c , (35) 

The slope m gives the critical exponent e and the Y axis intercept c is independent of /1. Fig. [6] shows the logarithm 
of Xq as a function of the logarithm of fi — hcep close to the CEP in QM model. Scaling is observed over several 



orders of magnitude. In Fig. 6(a) the [igep is approached from the lower /i side and we obtain a critical exponent 
0.6379 ± 0.0002 while the critical exponent e = m = 0.6648 ± 0.0001 in the result of Fig. |6(b)| whcn the 
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FIG. 6: (a) Shows the plot of the logarithm of \q as a function of the logarithm of // — hcep close to the CEP in QM model 
when the [icep is approached from the lower p side (b) shows the same plot as in Fig. a in QM model when the [icep is 
approached from the higher fi side. 
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FIG. 7: (a) Shows the plot of the logarithm of \q as a function of the logarithm of /j, — hcep close to the CEP in QMVT model 
when the ficEP is approached from the lower fi side (b) shows the same plot as in Fig.a in QMVT model when the (j,cep is 
approached from the higher \i side. 



IJ-CEP is approached from the higher /i side. The scaling starts around around log \fi — ^cep \ < — -5 in both the 
cases. These exponents show good agreement with the mean-field prediction e = 2/3. In Fig. 6(a)| the data fitted 
in a range 1.0 < log \(i — ^cep\ < 2.0 also shows a scaling kind of linear behaviour over one order of magnitude 
with a larger slope n — .6807 ± .0002 which changes to m — .6379 ± .0002 in the range —.5 < log \n — hcep\ < 1-0. 
When log\fi — Hcep\ ~ 2.0, we are very close to /i = on the temperature axis. The phase transition in the chiral 
limit at \x — in QM model is first order and it becomes crossover for the real life pion mass. In the quark mass 
(or the pion mass) and T plane at /i = 0, the first order transition line should change to crossover line through 
another second order critical end point as one increases the pion mass from zero to the experimental value. This 
linear behaviour in a range 1.0 < log|/i — /j.cep\ < 2.0 with larger slope may be due to the influence of another 
hidden CEP in the mass and temperature plane at \i = in the QM model. The critical exponent values obtained 
in PQM model calculation have been given in Table. 1, we can see that the presence of Polyakov loop potential in 
QM model, does not influence the value of critical exponents. Fig. [7] shows the plot of the logarithm of \q with 
respect to the logarithm of fj, — ficEP in the presence of fermionic vacuum fluctuations in QMVT model. In Fig. 7(a) 
when the hcep is approached from the lower /i side, we obtain a larger critical exponent e = m = 0.720 ± 0.00005 
in comparison to the corresponding result in the QM model. Due to the influence of fermionic vacuum fluctuation, 
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FIG. 8: (a) shows the plot of the logarithm of \q as a function of the logarithm of \i — hcep close to the CEP in PQMVT 
model when the \icep is approached from the lower \i side (b) shows the same plot as in Fig. a in PQMVT model when the 
HCep is approached from the higher /i side. 



Model 


H - fJ-CEP < 


M - ii-CEP > 


QM 
PQM 
QMVT 
PQMVT 


0.6379 ± 0.0002 
0.6309 ± 0.0001 
0.720 ± 0.00005 
0.725 ± 0.0002 


0.6648 ± 0.0001 
0.6668 ± 0.0001 
0.6938 ± 0.0002 
0.6886 ± 0.0004 



TABLE I: Critical exponents of the quark-number susceptibility in the QM,PQM,QMVT and PQMVT models for two different 
paths parallel to the chemical potential axis approaching the (icep from the lower /j, < [icep and higher fi > ficEpside. 



we find the presence of TCP in the chiral limit of QMVT model and it lies quite well within the R„ = 2 contour 



surrounding the CEP in the phase diagram in Fig. 1(a) This larger critical exponent may be the consequence of the 
modification of criticality around CEP due to the presence of TCP in its proximity. The scaling starts earlier when 
log [icep\ < 0.0 and we observe scaling over several orders of magnitude. This higher value of the critical exponent 
is close to the critical exponent calculated in Ref. [4j| where the effect of quantum fluctuations in the QM model were 
incorporated in the Proper-Time Renormalization Group (PTRG) approach. The critical exponents change in the 
range —0.5 < log \[i — hcep\ < -5 in Ref. [49j], from .77 for the one scaling regime starting after log — ^cep\ > -5 
to .74 for another scaling regime starting before log |/i — ^cep\ < — -5. Though we do not find analogous crossing 
behavior of the universality classes, the data points in our calculation show a bending trend and when we fit the 
data in a small range —0.5 < log |/i — /j,cep\ < -5, we find a higher slope e = m = 0.7664 ± 0.0002 as shown in 



the Fig. 7(a) In our calculation, the critical region of CEP is having a noticeable overlap with the critical region of 
TCP and this may be the reason of the bending trend in the data. If we approach the CEP from the higher /i side, 
we find smaller critical exponent e = m = 0.6938 ± 0.0002 in the result of Fig. |7(b)| In this case, the scaling starts 
around log |/i — ^cep\ < — .5 . It is pointed out that these exponents calculated in the presence of fermionic vacuum 
fluctuations in the QM model are different from the mean- field prediction e = 2/3. Similar results are found in Fig. [5] 
which shows the plot of the logarithm of \ q with respect to the logarithm of /i — [icep in the presence of fermionic 
vacuum fluctuations in PQMVT model. The presence of Polyakov loop compresses the width of critical region in 



PQMVT model but its effect on critical exponents is negligibly small as can be seen in Fig. [8(a)] and Fig. [8(b)] We 



point out that the recent calculation in 2+1 flavour model in Ref.[96[ does not report any modification of the mean 
field critical exponents under the infulence of fermionic vacuum term. The chiral crossover transition occurring at 
(i = on the temperature axis, will be faster and sharper due to the presence of s quarks in their 2+1 quark flavour 
calculation if we compare the present two quark flavour calculation with that of them. 

The critical exponents calculated in all the models are summarized and tabulated in Table 1. 
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FIG. 9: (a) Mass variations in the QM model are shown in this Fig. Two dotted lines in the rightmost part of the Fig. represent 
the temperature variations of m CT and at fi = 0, the dash dotted line in the middle part represents m CT while the line with 
small dash represents temperature variations at \jlcep = 151.7 and the solid line in the leftmost part of the Fig. represents 
m CT while the line with dash represents temperature variations at /i = 305 MeV. (b) Temperature variations of m a and 
in the influence of fermionic vacuum fluctuations, have been shown for the QMVT model calculations. Lines represent the 
same mass variations as in the left panel of the Fig. In QMVT model hcep = 299.35 MeV. 



In medium meson masses 



The critical fluctuations are also encoded in the variation of meson masses m w (T, fi) and m a (T, /i) as one passes 
through the chiral symmetry restoring phase transition. We will investigate and compare the 'in-medium' meson mass 
variations in QM, QMVT and PQM, PQMVT models in order to see the influence of fermionic vacuum fluctuations. 
The sigma and pion masses are calculated by determining the curvature of grand potential at the global minimum. 



dTTidTTi 



(36) 



ml(T,iJ,) 



dadc 



(37) 



The left panel of Fig. [9] shows the temperature variations of meson masses for fi = 0, fi = (J-cep and /i > fJ>cEP 
in QM model while the right panel shows the corresponding variations in QMVT model. In the chiral symmetry 
broken mesonic phase, the sigma mass always decreases with temperature. Sigma mass increases again at high 
temperatures signaling chiral symmetry restoration and it becomes degenerate with the increasing pion mass which 
does not vary much below the transition temperature. The degenerate meson masses increase linearly with T after 
the chiral symmetry restoration transition [49j has taken place. The temperature variations of m a and m w masses at 
fj, = in Fig. |9(bJ| are significantly modified due to the presence of fermionic vacuum fluctuations in QMVT model. 



If we compare these variations with the corresponding QM model temperature variations of masses in Fig. 9(a) 
find that the mass degeneration in m a and at = in Fig. |9(b)| becomes very smooth and it takes place at 
a higher temperature. Since the chiral crossover transition on the temperature axis at fi = is quite sharp and 
fast in QM model, the mass degeneration trend in m a and m T is also quite sharp and fast in Fig. 9(a) Fermionic 
vacuum fluctuations make the chiral crossover at /i = very smooth in QMVT model and this gets reflected also 
in the setting up of a very smooth mass degeneration trend at fi = in Fig. |9(b)| Long-wavelength fluctuations of 
the order parameter, characterize the second-order phase transitions. Since the chiral phase transition turns second 
order at the CEP, the sigma meson mass must vanish at the CEP because the effective potential completely flattens 
in the radial direction. Thus the sigma meson mass drops below the pion mass near the CEP. We know that the 
CEP in QMVT model gets shifted to a significantly higher chemical potential under the influence of fermionic vacuum 
fluctuation, hence the sigma mass becomes almost zero at ficEP — 299.35 MeV in QMVT model as shown in Fig. [9(b)] 
The sigma mass goes to zero only at y~cEP — 151.7 MeV in Fig. |9(b)| in QM model. The discontinuities in mass 
evolutions respectively in Fig. 9(a) and Fig. |9(b)| , signal a first-order phase transition at small temperatures when 
M > h^cep = 305 MeV in both the models QM as well as QMVT. 
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FIG. 10: (a) Mass variations in the PQM model are shown in this Fig. Two dotted lines in the rightmost part of the Fig. 
represent the temperature variations of m a and m-n at fi = 0, the dash dotted line in the middle part represents m a while the 
line with small dash represents m n temperature variations at hcep = 81.02 MeV and the solid line in the leftmost part of the 
Fig. represents m a while the line with dash represents m, temperature variations at fi — 305 MeV. (b) Temperature variations 
of m a and in the influence of fermionic vacuum fluctuations, have been shown for the PQMVT model calculations. Lines 
represent the same mass variations as in the left panel of the Fig. In PQMVT model hcep = 295.217 MeV. 



Fig. [10(a) 

TO, 



Finally we will be investigating how the fermionic vacuum term influences the emergence of masss degeneration 
trend in m a and as the chiral symmetry restoring transition takes place in the presence of Polyakov loop potential. 
Fig. 10(a) presents te mperatu re variations of meson masses at fi = 0, /x = [icep and /i > hcep in the in PQM model 
calculations while Fig. |10(b)| shows the corresponding mass variations in PQMVT model. Since the chiral crossover 
transition on the temperature axis at fj, = becomes very sharp and rapid in the PQM model due to the influence 
of Polyakov loop potential, the mass degeneration trend in m, a and at fj, = also becomes very sharp and fast in 
Here in the PQMVT model calculations also as in the case of QMVT model, the mass degeneration in 
and TOjr at /i = becomes very smooth in the influence of fermionic vacuum fluctuations as shown in Fig. |10(b)| 
Further this mass degeneration takes place at a temperature which is highest of the corresponding mass degeneration 
temperatures in other models. This happens because, the chiral crossover transition at /i = occurs at a temperature 
which is highest in the PQMVT model calculations, due to the combined effect of Polyakov loop potential and 
fermionic vacuum term. We point out that the sigma meson mass neither vanishes nor becomes very close to zero 
at CEP as shown in Fig. 10(a) and Fig. |10(bJ| respectively for PQM and PQMVT model temperature variations. In 
PQM model, m a in its temperature variation reaches the minimum value of 128 MeV as shown by the dash dotted 
line in Fig. 



io( a ; 



while in PQMVT model the minimum value reached in the temperature variation of m a is 76.0 MeV 
as evident from the dash dotted line in Fig. |10(b)| It means that the effective potential of PQM and PQMVT model 
does not completely flattens in the radial direction at the CEP. The Polyakov loop expectation value is a scalar field 
which mixes up with the chiral order parameter and this effect hampers the flattening of the PQM model effective 
potential in the radial direction (i.e. the direction of a field) at the CEP. In PQMVT model, this effect seems to be 
considerably remedied by the presence of the fermionic vacuum term and the minimum value of m a becomes 76.0 
MeV in its temperature variation. Here also, the discontinuities in mass evolutions respectively in Fig. 10(a) and 
Fig. 1 10 (b) | , signal a first-order phase transition at small temperatures when \x > [icep — 305 MeV in both the models 
PQM as well as PQMVT. 



IV. SUMMARY AND CONCLUSION 



In the beginning of the present work, we computed the phase diagrams and pin-pointed the CEP positions in the 
H and T plane of all the four two quark flavor models QMVT, PQMVT, QM and PQM for the real life explicit chiral 
symmetry breaking with the experimental value of pion mass. We obtained the phase diagrams with zero pion mass 
also for the chiral limit in the QMVT and PQMVT models and located the respective positions of TCP. Since the 
presence or absence of TCP in the phase diagram and its distance from CEP in a model calculation, influences the 
nature of critical fluctuations around CEP, we quantified the proximity of TCP to the CEP in the phase diagrams 
obtained for QMVT and PQMVT model calculations. 

The QMVT model CEP is positioned at ^ CjBP =299.35 MeV and T CBP =32.24 MeV and it shifts to higher tern- 
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perature Tcep=83.0 MeV and lower chemical potential /j.cep—295 .217 MeV in PQMVT model due to the effect of 
Polyakov loop. In QM model, the CEP is located at T CBP =102.09 MeV and fj, CEP =151.7 MeV and again in the 
influence of Polyakov loop potential, it shows considerable shift towards the temperature axis in PQM model and gets 
located at T C£ ;p=166.88 MeV and ^ CEP =81.02 MeV. If we compare the location of CEP in QM and PQM models to 
the location of CEP in QMVT and PQMVT models, we find a considerably significant shift of CEP to large chemical 
potential and small temperature values for QMVT and PQMVT models due to the robust influence of fermionic 
vacuum term presence in the effective potential. In chiral limit of PQMVT model, the tricritical point(TCP) exists at 
T t =137.09 MeV and ^=240.14 MeV and the TCP in QMVT model is found at T t =69.06 MeV and Mt =263.0 MeV. 
The proximity of TCP to the CEP has been quantified by plotting the constant normalized quark-number susceptibil- 
ity (R q =2) contours around CEP in the phase diagrams of PQMVT and QMVT models. The second cumulant of the 
net quark number fluctuations on these contours is double to that of the free quark gas value and such enhancements 
are the signatures of CEP for the heavy-ion collision experiments. The TCP location is quite well inside the R q =2 
contour on the phase diagram of both the models QMVT as well as PQMVT. 

The CEP in phase diagrams is pin-pointed by tracking down the divergence in the quark number susceptibilities 
and scalar susceptibilities which show significant enhancement in a region around the CEP in the [i and T plane 
in comparison to their respective values for the free quark gas. In order to determine the shape of the critical 
region around CEP, we have plotted different contour regions having different constant values of properly normalized 
quark number susceptibility ratio (R q ) and scalar susceptibility ratio (Rs). The different shapes of these contours as 
calculated in various models, throw light on the nature of criticality around CEP in those models. We have plotted 
the contours with three different values of quark number susceptibility ratio R q = 2,3 and 5, in the /i and T plane 
relative to the CEP. If we compare the contours obtained in the PQM model to the contours in pure QM model, we 
conclude that the presence of Polyakov loop potential, compresses the critical region particularly in the T direction. 
Since the chiral crossover transition becomes faster and sharper due to the presence of Polyakov loop potential, the 
critical region in the T direction gets significantly compressed. The analysis of the shape of R q = 2, 3 and 5 contours 
in the QMVT and PQMVT models, tells us that the size of critical region is increased in a direction perpendicular to 
the crossover line due to the influence of the fermionic vacuum fluctuations. This effect is less pronounced in PQMVT 
model because of the compression of critical region width due to the presence of Polyakov loop potential while the 
QMVT model contours show a robust increase in the width of the critical region in perpendicular direction to the 
crossover transition line . However, the extent and size of critical region in the PQMVT model is noticeably larger in 
both the directions fi as well as T compared to that of QMVT model results. In the presence of fermionic vacuum term, 
CEP gets located at larger chemical potentials in QMVT/PQMVT models. Since the quark determinant gets modified 
mostly at moderate chemical potentials by the presence of Polyakov loop potential and further in its influence, the 
PQMVT model CEP shifts to a higher critical temperature when compared to the CEP in QMVT model, we obtain 
an enhancement of the critical region in PQMVT model. 

We have plotted three contours around the CEP in PQM and QM models also for the normalized scalar susceptibility 
ratios R s = 10, 15 and 25. The shape of R s = 10 contour in PQM model is compressed in comparison to the pure 
QM model contours. The R s = 25 contour is not obtained in PQM model calculation because the minimum value 
of a meson mass does not fall below 100 MeV. Since the value of m a falls very rapidly and sharply from 500 MeV 
to 128 MeV, we obtain a very thin and small contour region even for R s — 15. We get all the contours regions 
for R s = 10, 15 and 25 with well defined size in the QM model calculations because the m a variation is smoother 
and slower in comparison to the corresponding PQM model results and further the minimum in the vn a variation 
approaches almost zero value in QM model. Further we do not find the closure of R s = 10 contour on the temperature 
axis at fi = MeV in both the models QM as well as PQM. We obtain quite well defined and closed contour regions 
for R s = 10, 15 and 25 in QMVT and PQMVT model calculations which again become broader in the direction 
perpendicular to the crossover line due to the presence of fermionic vacuum fluctuations. For scalar susceptibility 
also, the critical region gets elongated in the phase diagram and x<r is enhanced in the direction parallel to the 
first-order transition line. Here also the presence of Polyakov loop potential in the PQMVT model, leads to the 
compression in the width of critical region around CEP. The fermionic vacuum fluctuations, make the chiral crossover 
transition very smooth while the Polyakov loop potential makes it sharper and faster and these opposite effects give 
a typical shape to the quark number susceptibility contours in the PQMVT model. Similar effects can be seen in the 
scalar susceptibility contours also. In the influence of fermionic vacuum fluctuations only, the x<r contours in the pure 
QMVT model are broader and rounded. 

In order to further investigate the nature of criticality in two flavour calculations, we have numerically evaluated 
the critical exponents of the quark number susceptibility Xq m QM, PQM, QMVT and PQMVT models. In these 
investigations, the critical hcep at fixed critical temperature Tqep is approached from the lower as well as higher \x 
sides in a path parallel to the /x-axis in the (T, /i)-plane. The calculation of the critical exponents has been done using 
the linear logarithmic fit. If the \icep is approached in QM model from the lower /j, side, we obtain a critical exponent 
equal to e = 0.6379 ± 0.0002 while the critical exponent is e = 0.6648 ± 0.0001 when the \xcep is approached from the 
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higher /i side. The scaling starts around around log \ji — fJ>CEp\ < — -5 in both the cases. These exponents show good 
agreement with the mean-field prediction e = 2/3. The influence of Polyakov loop potential on the calculated values 
of critical exponents in PQM model, is negligible and we obtain similar critical exponents as evaluated in QM model. 

If the hcep hi QMVT model calculation is approached from the lower side, we obtain a critical exponent 
e = m = 0.720 ± 0.00005 which is larger in comparison of the corresponding critical exponent e = 0.6379 ± 0.0002 
evaluated in the QM model. Due to the influence of fermionic vacuum fluctuation, we find the presence of TCP in 
the chiral limit of QMVT model and it lies quite well within the R q — 2 contour surrounding the CEP in the phase 
diagram. This larger critical exponent may be the consequence of the modification of criticality around CEP due to 
the presence of TCP in its proximity. The scaling starts earlier when log — ^cep\ < 0.0 and we observe scaling over 
several orders of magnitude. We obtain smaller critical exponent e = 0.6938 ± 0.0001 when the ^cep is approached 
from the higher [i side in PQMVT model. It is pointed out that these exponents calculated in the presence of fermionic 
vacuum fluctuations in the QM model are different from the mean-field prediction e = 2/3. The presence of Polyakov 
loop compresses the width of critical region in PQMVT model but its effect on critical exponents is negligible. The 
critical exponents obtained in PQM model calculations are also similar to the critical exponents of QM model. 

Since the critical fluctuations are encoded in the variation of meson masses (T, \i) and ra a (T, fi) as one passes 
through the chiral symmetry restoring phase transition, we have also investigated and compared the 'in-medium' 
meson mass temperature variations for \i = 0, yU = [icep and > [icep in QM, QMVT and PQM, PQMVT model 
calculations. If we compare the temperature variations of masses at = in QM and QMVT model calculations, we 
find that the mass degeneration in m CT and at [i — becomes very smooth in QMVT model and it takes place at 
a higher temperature. Since the sharper and faster chiral crossover transition occurring on the temperature axis at 
/i = in QM model, becomes very smooth in QMVT model under the influence of fermionic vacuum fluctuations, the 
sharper mass degeneration trend in m a and in QM model also becomes a very smooth mass degeneration trend in 
QMVT model. The sigma meson mass must vanish at the CEP since the chiral phase transition turns second order 
at this point and the effective potential completely flattens in the radial direction. In our QM and QMVT model 
calculations, we have shown that the sigma meson mass becomes almost zero at /i = (J-cep- in both the models. It has 
also been shown that the discontinuities in mass evolutions, signal a first-order phase transition at small temperatures 
when \x > \x C ep = 305 MeV in both the models QM as well as QMVT. 

Finally our investigation gets concluded by the study of the influence of Polyakov loop potential on the emergence 
of masss degeneration trend in m a and m, as the chiral symmetry restoration takes place in the presence of fermionic 
vacuum fluctuation. Since the chiral crossover transition on the temperature axis at [i = becomes very sharp and 
rapid in the PQM model due to the influence of Polyakov loop potential, the mass degeneration trend in m a and 
at,u = also becomes very sharp and fast in PQM model. Here in the PQMVT model calculations also as in the case 
of QMVT model, the mass degeneration in m CT and at yu = becomes very smooth in the influence of fermionic 
vacuum fluctuations. Further this mass degeneration happens at a temperature which is highest of the corresponding 
mass degeneration temperatures in other models. In the PQMVT model calculations, we obtain highest temperature 
for chiral crossover transition on the temperature axis at fj, = 0. The combined effect of Polyakov loop potential and 
fermionic vacuum term is responsible for this. We point out that the sigma meson mass neither vanishes nor becomes 
very close to zero at [i — ^cep in PQM and PQMVT model temperature variations. In PQM model, m a in its 
temperature variation reaches the minimum value of 128.0 MeV while the minimum value reached in the temperature 
variation of m a in PQMVT model is 76.0 MeV. It means that the effective potential of PQM and PQMVT model 
does not completely flattens in the radial direction at the CEP. The Polyakov loop expectation value is a scalar field 
which mixes up with the chiral order parameter and this effect hampers the flattening of the PQM model effective 
potential in the radial direction (i.e. the direction of afield) at the CEP. In the PQMVT model, this effect seems 
to be considerably remedied by the presence the fermionic vacuum term and we get minimum in m a temperature 
variation at 76.0 MeV. We have also shown the discontinuities in mass evolutions which signal a first-order phase 
transition at small temperatures when [i > [icep = 305 MeV in both the models PQM as well as PQMVT. 
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